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The  following  three  reports  are  being  distributed  together  and  reprints 
of  the  papers  previously  published  with  their  texts,  in  scientific  journals, 
are  used  for  the  purpose.  Circumstantial  reasons,  primarly  cost,  made  this 
procedure  advisable.  One  of  the  reports  (56)  deals  with  the  optimization 
of  the  shape  of  an  origionally  circular  hole  present  in  a  square  plate, 
loaded  uniformly  and  uniaxially  in  its  plane.  A  second  report  (57)  deals 
with  the  optimization  of  the  shape  of  the  hole,  when  the  hole  is  present 
in  a  cylindrical  shell  subjected  to  axial  tension.  The  third  report  (55) 
deals  with  the  necessary  background  for  the  optimization  of  the  hole  in 
the  square  plate. 

The  abstracts  of  the  previously  published  reports,  prepared  with  O.N.R. 
support,  are  presented  only  once,  before  the  text  of  the  three  present  reports. 
The  O.N.R.  reports  distribution  list  is  also  presented  only  once,  at  the 
end.  A  Report  Documentation  Page  follows  each  individual  report. 
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Previous  Technical  Reports  to  the  Office  of  Naval  Research 


1.  A.  J.  Durelli,  "Development  of  Experimental  Stress  Analysis  Methods  to 
Determine  Stresses  and  Strains  in  Solid  Propellant  Grains" --June  1962. 
Developments  in  the  manufacturing  of  grain-propellant  models  are 
reported.  Two  methods  are  given:  a)  cementing  routed  layers  and 

b)  casting. 

2.  A.  J.  Dxirelli  and  V.  J.  Parks,  "New  Method  to  Determine  Restrained 
Shrinkage  Stresses  in  Propellant  Grain  Models" — October  1962. 

The  birefringence  exhibited  in  the  curing  process  of  a  partially 
restrained  polyurethane  rubber  is  used  to  determine  the  stress  associated 
with  restrained  shrinicage  in  models  of  solid  propellant  grains  partially 
bonded  to  the  case. 

3.  A.  J.  Durelli,  "Recent  Advances  in  the  Application  of  Photoelasticity  in 
the  Missile  Industry" — October  1962. 

Two-  and  three-dimensional  photoelastic  analysis  of  grains  loaded  by 
pressure  and  by  temperature  are  presented.  Seme  applications  to  the 
optimization  of  fillet  contours  and  to  the  redesign  of  case  joints  are 
also  included. 

4.  A.  J.  Durelli  and  V.  J.  Parks,  "Experimental  Solution  of  Some  Mixed 
Boundary  Value  Problems"— April  1964. 

Means  of  applying  known  displacements  and  known  stresses  to  the  boundaries 
of  models  used  in  experimental  stress  analysis  are  given.  The  applica¬ 
tion  of  some  of  these  methods  to  the  analysis  of  stresses  in  the  field 
of  solid  propellant  grains  is  illustra-ced.  The  presence  of  the  "pinching 
effect"  is  discussed. 

5.  A.  J.  Durelli,  "Brief  Review  of  the  State  of  the  Art  and  Expected  Advance 
ip  Experimental  Stress  and  Strain  Analysis  of  Solid  Propellant  Grains" — 
April  1964. 

A  brief  review  is  made  of  the  state  of  the  experimental  stress  and  strain 
analysis  of  solid  propellant  grains.  A  discussion  of  the  prospects  for 
the  next  fifteen  years  is  added. 

t 

6.  A.  J.  Durelli,  "Experimental  Strain  and  Stress  Analysis  of  Solid  Propellant 
Rocket  Motors" — March  1965. 

A  review  is  made  of  the  experimental  methods  used  to  strain-analyze  solid 
propellant  rocket  motor  shells  and  grains  when  subjected  to  different 
loading  conditions.  Methods  directed  at  the  determination  of  strains  in 
actual  rockets  are  included. 

7.  L.  Ferrer,  V.  J.  Parks  and  A.  J.  Durelli,  "An  Experimental  Method  to  Analyze 
Gravitational  Stresses  in  Two-Dimensional  Problems"— October  1965. 
Photoelasticity  and  moire  methods  are  used  to  solve  two-dimensional  problems 
in  which  gravity-stresses  are  present. 
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8. 


A.  J.  Durelii,  V.  J.  Parks  and  C.  J.  del  Rio,  "Stresses  in  a  Square  Slab 
Bonded  on  One  face  :o  a  Rigid  Plate  and  Shrunk"— November  1965. 

A  square  epoxy  slab  was  bonded  to  a  rigid  plate  on  one  of  its  faces  in 
the  process  of  curing.  In  the  same  process  the  photoelastic  effects 
associated  with  a  state  of  restrained  shrinkage  were  "frozen-in." 
Three-dimensional  photoelasticity  was  used  in  the  analysis. 

9.  A.  J.  Durelli,  V.  J.  Parks  and  C.  J.  del  Rio,  "Experimental  Determination 
of  Stresses  and  Displacements  in  Thick-Wall  Cylinders  of  Complicated 
Shape" — April  1966. 

Photoelasticity  and  moire  are  used  to  analyze  a  three-dimensional  rocket 
shape  with  a  star  shaped  core  subjected  to  internal  pressure. 

10.  V.  J.  Parks,  A.  J.  Durelli  and  L.  Ferrer,  "Gravitational  Stresses 
Determined  Using  Immersion  Techniques"— July  1966. 

The  methods  presented  in  Technical  Report  No.  7  above  are  extended  to 
three-dimensions.  Immersion  is  used  to  increase  response. 

11.  A.  J.  Durelli  and  V.  J.  Parks.  "Experimental  Stress  Analysis  of  Loaded 
Boundaries  in  Two-Dimensional  Second  Boundary  Value  Problems" — 

February  1967. 

The  pinching  effect  that  occurs  in  two-dimensional  bonding  problems, 
noted  in  Reports  2  and  4  above,  is  analyzed  in  some  detail. 

12.  A.  J.  Durelli,  V.  J.  Parks,  H.  C.  Feng  and  F.  Chiang,  "Strains  and 
Stresses  in  Matrices  with  Inserts," —  May  1967. 

Stresses  and  strains  along  the  interfaces,  and  near  the  fiber  ends,  for 
different  fiber  end  configiirations,  are  studied  in  detail. 

13.  A.  J.  Durelli,  V.  J.  Parks  and  S.  Uribe,  "Optimization  of  a  Slot  End 
Configuration  in  a  Finite  Plate  Subjected  to  Uniformly.  Distributed 
Load,"— June  1967. 

Two-dimensional  photoelasticity  was  used  to  study  various  elliptical  ends 
to  a  slot,  and  determine  which  would  give  the  lowest  stress  concentration 
for  a  load  normal  to  the  slot  length. 

14.  A.  J.  Durelli,  V.  J.  Parks  and  Han-Chow  Lee,  "Stresses  in  a  Split 
Cylinder  Bonded  to  a  Case  and  Subjected  to  Restrained  Shrinkage," — 

Janiiary  1968. 

A  three-dimensional  photoelastic  study  that  describes  a  method  and 
shows  results  for  the  stresses  on  the  free  boundaries  and  at  the 
bonded  interface  of  a  solid  propellant  rocket. 

15.  A.  J.  Durelli,  "Experimental  Stress  Analysis  Activities  in  Selected 
European  Laboratories" — August  1968. 

This  report  has  been  written  following  a  trip  conducted  by  the  author 
through  several  European  countries.  A  list  is  given  of  many  of  the 
laboratories  doing  important  experimental  stress  analysis  work  and  of 
the  people  interested  in  this  kind  of  work.  An  attempt  has  been  made 
to  abstract  the  main  characteristics  of  the  methods  used  in  some  of 
the  countries  visited. 
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V.  J.  Parks,  A.  J.  Dureili  and  L.  Ferrer,  "Constant  Acceleration  Stresses 
in  a  Composite  Body" — October  1968. 

Use  of  the  immersion  analogy  to  determine  gravitational  stresses  in 
two-drnensional  bodies  made  of  materials  with  different  properties. 

,  17.  A.  J.  Dureili,  J.  A.  Clark  and  A.  Kochev,  "Experimental  Analysis  of  High 

Frequency  Stress  Waves  in  a  Ring"— October  1968. 

A  method  for  the  complete  experimental  determination  of  dynamic  stress 
distributions  in  a  ring  is  demonstrated.  Photoelastic  data  is  supple¬ 
mented  by  measurements  with  a  capacitance  gage  used  as  a  dynamic  lateral 
extensometer . 

18.  J.  A.  Clark  and  A.  J.  Dureili,  "A  Modified  Method  of  Holographic  Inter¬ 
ferometry  for  Static  and  Dynamic  Photoelasticity" — April  1968. 

A  simplified  absolute  retardation  approach  to  photoelastic  analysis  is 
described.  Dynamic  isopachics  are  presented. 

19.  J.  A.  Clark  and  A.  J.  Dureili,  "Photoelastic  Analysis  of  Flexural  Waves 
in  a  Bar"— May  1969. 

A  complete  direct,  full-field  optical  determination  of  dynamic  stress 
distribution  is  illustrated.  The  method  is  applied  to  the  study  of 
flexural  waves  propagating  in  a  urethane  rubber  bar.  Results  are 
compared  with  approximate  theories  of  fle.xural  waves. 

20.  J.  A.  Clark  and  A.  J.  Dureili,  "Optical  Analysis  of  Vibrations  in 
Continuous  Media"— June  1969. 

Optical  methods  of  vibration  analysis  are  described  which  are  independent 
of  assumptions  associated  with  theories  of  wave  propagation.  Methods  are 
illustrated  with  studies  of  transverse  waves  in  prestressed  bars,  snap 
loading  of  bars  and  motion  of  a  fluid  surrounding  a  vibrating  bar. 

21.  V.  J.  Parks,  A.  J.  Dureili,  K.  Chandrashekhara  and  7.  L.  Chen,  "Stress 
Distribution  Around  a  Circular  Bar,  with  Flat  and  Spherical  Ends, 

Embedded  in  a  Matrix  in  a  Triaxial  Stress  Field"— July  1969. 

A  Three-dimensional  photoelastic  method  to  determine  stresses  in  composite 
materials  is  applied  to  this  basic  shape.  The  analyses  of  models  with 
differe.nt  loads  are  combined  to  obtain  stresses  for  the  triaxial  cases. 

22.  A.  J.  Dureili,  V.  J.  Parks  and  L.  Ferrer,  "Stresses  in  Solid  and  Hollcw 
Spheres  Subjected  to  Gravity  or  to  Normal  Surface  Tractions"— 

October  1969. 

The  method  described  in  Report  No.  10  above  is  applied  to  two  specific 
problems.  An  approach  is  suggested  to  extend  the  solutions  to  a  class 
of  surface  traction  problems. 

23.  J.  A.  Clark  and  A.  J.  Dureili,  "Separation  of  Additive  and  Subtractive 
Moir4  Patterns"— December  1969. 

A  spatial  filtering  technique  for  adding  and  subtracting  images  of  several 
gratings  is  described  and  employed  to  determine  the  whole  field  of 
Cartesian  shears  and  rigid  rotations. 
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2“.  S.  Sanford  and  A.  J.  Durelii,  ''Interpretation  of  Fringes  in  Stress- 
Holo- Interf erometry" --T -iy  137  7 . 

Errors  associated  with  interpreting  stress-hoio-interferonetry  patterns 
as  the  superposition  of  isopachics  (with  half  order  fringe  shifts)  and 
isochromatics  are  analysed  theoretically  and  illustrated  with  computer 
generated  holographic  interference  patterns. 

25.  J.  A.  Clark,  A.  J.  Durelii  and  ?.  A.  Laura,  "On  the  Effect  of  Initial 
Stress  on  the  Propagation  of  Flexural  Waves  in  Elastic  Rectangular 
3ars" — December  1370. 

Experimental  analysis  of  the  propagation  of  flexural  waves  in  prismatic, 
elastic  bars  with  and  without  prestressing.  Tne  effects  of  prestressing 
by  axial  tension,  axial  compression  and  pure  bending  are  illustrated. 

2o.  A.  J.  Durelii  and  J.  A.  Clark,  "Experimental  Analysis  of  Stresses  in  a 
Buoy-Cable  System  Using  a  Birefringent  Fluid" — February  1971. 

An  extension  of  the  method  .of  photoviscous  analysis  is  presented  which 
permits  quantitative  studies  of  strains  associated  with  steady  state 
vibrations  of  immersed  structures.  The  method  is  applied  in  an 
investigation  of  one  form  of  behavior  of  buoy-cable  systems  loaded  by 
the  action  of  surface  waves. 

27.  A.  J,  Durelii  and  T.  L.  Chen,  "Displacements  and  Finite-Strain  Fields  in 
a  Sphere  Subjected  to  Large  Defortr.ations"--February  1972. 

Displacements  and  strains  (ranging  from  0.001  to  0.50)  are  determined  in 

a  polyurethane  sphere  subjected  to  several  levels  of  diametral  compression. 
A  500  lines -per -inch  grating  was  embedded  in  a  meridian  plane  of  the 
sphere  and  moire  effect  produced  with  a  non-deformed  master.  The  maximum 
applied  vertical  displacement  reduced  the  diameter  of  the  sphere  by  27 
per  cent. 

28.  A.  J.  Durelii  and  S.  Machida,  "stresses  and  Strain  in  a  Disk  with  Variable 
Modulus  of  Elasticity" — March  1972 

A  transparent  material  with  variable  modulus  of  elasticity  has  been 
manufactured  chat  exhibits  good  pbocoelascic  properties  and  can  also  be 
strain  analyzed  by  moir^.  The  results  obtained  suggests  Chat  the  stress 
distribution  in  the  disk  of  variable  E  i.s  practically  Che  same  as  Che 
stress  distribution  in  the  homogeneous  disk.  It  also  indicates  chat  Che 
strain  fields  in  both  cases  are  very  different,  but  that  It  is  possible, 
approximately,  to  obtain  Che  stress  field  from  the  strain  field  using  Che 
value  of  E  at  every  point,  and  Hooke's  law. 

29.  A.  J.  Durelii  and  .1.  Ruitrago,  "State  of  Stress  and  Strain  in  a  Rectangular 
Belt  Hulled  Over  a  Cylindrical  Pulley" — June  1972. 

Two-  and  three-dimensional  photoelascicitv  as  well  as  electrical  strain 
gages,  dial  gages  and  micrometers  are  used  to  determine  the  stress  distri¬ 
bution  in  a  belt-pulley  system.  Contact  .and  tangential  stress  for  various 
contact  angles  and  friction  coefficients  are  given. 
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30.  T.  L.  Chen  and  A.  J.  Ourelli,  "Stress  Field  in  a  Sphere  Subjected  to 
large  Deformations"— June  1972. 

Strain  fields  obtained  in  a  sphere  subjected  to  large  diametral  compressions 
from  a  previous  paper  were  converted  into  stress  fields  using  two  approaches. 
First,  the  concept  of  strain-energy  function  for  an  isotropic  elastic 
body  was  used.  Then  the  stress  field  was  determined  with  the  Hookean 
type  natural  stress-natural  strain  relation.  The  results  so  obtained 
were  also  compared. 

31.  A.  J.  Ourelli,  V.  J.  Parks  and  H.  M.  Kasseem,  "Helices  Under  Load" — 

July  1973. 

Previous  solutions  for  the  case  of  close  coiled  helical  springs  and  for 
helices  made  of  thin  bars  are  extended.  The  complete  solution  is 
presented  in  graphs  for  the  use  of  designers.  The  theoretical  development 
is  correlated  with  expe-  ^nts. 

32.  T.  L.  Chen  and  A.  J.  Ourelli,  "Displacements  and  Finite  Strain  Fields  in 
a  Hollow  Sphere  Subjected  to  Large  Elastic  Deformations"— September  1973. 

The  same  methods  described  in  No.  27,  were  applied  to  a  hollow  sphere 
with  an  inner  diameter  one  half  the  outer  diameter.  The  hollow  sphere 
was  loaded  up  to  a  strain  of  30  per  cent  on  the  meridian  plane  and  a 
reduction  of  the  diameter  by  20  per  cent. 

33.  A.  J.  Durelli,  H.  H.  Hasseem  and  V.  J.  Parks,  "New  Experimental  Method 
in  Three-Dimensional  Elastostatics" — December  1973. 

A  new  material  is  reported  which  is  unique  among  three-dimensional 
stress-freezing  materials,  in  that,  in  its  heated  (or  rubbery)  state 
it  has  a  Poisson's  ratio  which  is  appreciably  lower  than  0.5.  For  a 
loaded  model,  made  of  this  material,  the  unique  property  allows  the 
direct  determination  of  stresses  from  strain  measurements  ’taken  at 
interior  points  in  the  model. 

34.  J.  Wolak  and  V.  J.  Parks,  "Evaluation  of  Large  Strains  in  Industrial 
Applications" — April  1974. 

It  was  shown  that  Mohr's  circle  permits  the  transformation  of  strain  from 
one  axis  of  reference  to  another,  irrespective  of  the  magnitude  of  the 
strain,  and  leads  to  the  evaluation  of  the  principal  strain  components 
from  the  measurement  of  direct  strain  in  three  directions. 

35.  A.  J.  Durelli,  "Experimental  Stress  Analysis  Activities  in  Selected 
European  Laboratories"— April  1975. 

Continuation  of  Report  No.  15  after  a  visit  to  Belgium,  Holland,  Germany, 

France,  Turkey,  England  and  Scotland. 

36.  A.  J.  Durelli,  V.  J.  Parks  and  J.  0.  Biihler-Vidal,  "Linear  and  Non-linear 
Elastic  and  Plastic  Strains  in  a  Plate  with  a  Big  Hole  Loaded  Axially  in 
its  Plane"— July  1975. 

Strain  analysis  of  the  ligament  of  a  plate  with  a  big  hole  indicates.«that 
both  geome-oric  and  material  non-linearity  may  take  place.  The  strain 
concentration  factor  was  found  to  vary  from  1  to  2  depending  on  the  level 
of  deformation. 
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37.  A.  J.  Dur«lii,  V.  Paviin,  J.  0.  Biihier-Vidai  and  G.  Ome,  "Elastostatics 
of  a  Cubic  Box  Sab j  acted  to  Concentrated  Loads  "»August  1975. 

Analysis  of  experimental  strain,  stress  and  deflection  of  a  cubic  box 
subjected  to  concentrated  loads  applied  at  the  center  of  two  opposite 
faces.  The  ratio  between  the  inside  span  and  the  wall  thickness  was 
varied  between  approximately  5  and  121. 

38.  A.  J.  Iha:*elii,  V.  J.  Parks  and  J.  0.  BUhler-Vidal,  "Elastostatics  of 
Cubic  Boxes  Subjected  to  Pressure" — March  1976. 

Experimental  analysis  of  strain,  stress  and  deflections  in  a  cubic  box 
subjected  to  either  internal  or  external  pressure.  Inside  span-to-wall 
thickness  ratio  varied  from  5  to  14. 

39.  Y.  Y.  Hung,  J.  D.  Hovanesian  and  A.  J.  Durelli,  "New  Optical  Method  to 
Determine  Vibration-Induced  Strains  with  Variable  Sensitivity  After 
Recoirding"— November  1976. 

A  steady  state  vibrating  object  is  illuminated  with  coherent  light  and 
its  image  slightly  misfocused.  The  resulting  specklegrara  is  "time- 
integrated"  as  when  Fouri^  filtered  gives  derivatives  of  the  vibrational 
amplitude. 

40.  Y.  Y.  Hung,  C.  Y.  Liang,  J.  D.  Hovanesian  and  A.  J.  Durelli,  "Cyclic 
Stress  Studies  by  Time-Averaged  Photoelasticity"— November  1976. 
"Time-averaged  isochromatics"  are  formed  when  the  photographic  film  is 
exposed  for  more  than  one  period.  Fringes  represent  amplitudes  of  the 
oscillating  stress  according  to  the  zeroth  order  Bessel  function. 

41.  Y.  Y.  Hung,  C.  Y.  Liang,  J.  D.  Hovanesian  and  A.  J.  Durelli,  "Time- 
Averaged  Shadow  Moir^  Method  for  Studying  Vibrations"— November  1976. 
Time-averaged  shadow  raoir€  petmiits  the  determination  of  the  amplitude 
distribution  of  the  deflection  of  a  steady  vibrating  plate-. 

42.  J.  Buitrago  and  A.  J.  Durelli,  "On  the  Interpretation  of  Shadow-Mo ird 
Fringes"— April  1977. 

Possible  rotations  and  translations  of  the  grating  are  considered  in  a 
general  expression  to  interpret  shadow-moir^  fringes  and  on  the 
sensitivity  of  the  method.  Application  to  an  inverted  perforated  tube. 

43.  J.  der  Hovanesian,  "18th  Polish  Solid  Mechanics  Conference."  Published  in 
European  Scientific  Notes  of  the  Office  of  Naval  Research,  in  London, 
England,  Dee.  31,  1976. 

Comments  on  the  planning  and  organization  of,  and  scientific  content  of 
paper  presented  at  the  18th  Polish  Solid  Mechanics  Conference  held  in 
Hisla-Jawornik  from  September  7-14,  1976. 

44.  A.  J.  Durelli,  "The  Difficult  Choice," — May  1977, 

The  advantages  and  limitations  of  methods  available  for  the  analyses 
of  displacements,  strain,  and  stresses  are  considered.  Comments  are 
made  on  several  theoretical  approaches,  in  particular  approximate 
methods,  and  attention  is  concentrated  on  experimental  methods:  photo¬ 
elasticity,  moire,  brittle  and  photoelastic  coatings,  gages,  grids, 
holography  and  speckle  to  solve  two-  and  three-dimensional  problems  in 
elasticity,  plasticity,  dynamics  and  anisotropy. 
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45>  C.  Y.  Liang,  Y.  Y.  Hung,  A.  J.  Ourelli  and  J.  0.  Hovaneslan, 

"Direct  Oecermlnacion  of  Flexural  Strains  in  Places  Using  Projected 
Gratings ,  "—June  1977. 

The  oechod  requires  the  rotation  of  one  photograph  of  the  deformed 
gracing  over  a  copy  of  itself.  The  moird  produced  yields  strains  by 
optical  double  differentiation  of  deflections.  Applied  to  projected 
gratings  Che  idea  permits  the  study  of  plates  subjected  to  much  larger 
deflections  chan  the  ones  chat  can  be  studied  with  holograms. 


46.  A.  J.  Durelli,  K.  Brown  and  P.  Yee,  "Optimization  of  Geometric 
Ois continuities  in  Stress  Fields" — March  1978. 

The  concept  of  "coefficient  of  efficiency"  is  introduced  to  evaluate 
Che  degree  of  optimization.  An  ideal  design  of  the  inside  boundary  of 
a  cube  subjected  to  diametral  coapression  is  developed  which  decreases 
its  maximum  stress  by  25Z,  at  the  time  it  also  decreases  its  weight  by 
lOZ.  The  efficiency  coefficient  is  Increased  from  0.59  to  0.95. 

Tests  with  a  brittle  material  show  an  increase  in  strength  of  20Z.  An 
ideal  design  of  the  boundary  of  the  hole  in  a  place  subjected  to  axial 
load  reduces  the  maximum  stresses  by  26Z  and  increases  the  coefficient 
of  efficiency  from  0.54  to  0.90. 

47.  J.  D.  Hovaneslan,  Y.  Y.  Hung  and  A.  J.  Durelli,  "New  Optical  Method 

to  Determine  Vibration-Induced  Strains  With  Variable  Sensitivity  After 
Recording" — May  1978. 

A  steady-state  vibrating  object  is  illuminated  with  coherent  light  and 
its  image  is  slightly  mis  focused  in  the  film  plane  of  a  camera.  The 
resulting  processed  film  is  called  a  "time-integrated  specklegram. " 

When  the  specklegram  is  Fourier ‘filtered,  it  e.xhibics  fringes  depicting 
derivatives  of  the  vibrational  amplitude.  The  direction  of  the  spatial 
derivative,  as  well  as  the  fringe  sensitivity  may  be  easily  and  continu¬ 
ously  varied  during  the  Fourier  filtering  process.  This  new  method  is 
also  much  less  demanding  than  holographic  interferometry  with  respect  to 
vibration  isolation,  optical  set-up  cine,  illuminating  source  coherence, 
required  film  resolution,  etc. 


48.  Y.  Y.  Ilunv  A.  .1.  Durelli,  imulcaneous  Determination  of  Three 
Strain  Components  in  Speckle  Interferometry  Uslnp  a  Multiple  Image 
She.'irinp  C.imera," — September  1978 

This  paper  describes  a  multiple  image-shearing  camera.  Incorporating 
coherent  light  illumination,  Che  camera  serves  as  a  multiple  shearing 
speckle  interferometer  which  measures  the  derivatives  of  surface 
displacements  with  respect  to  three  directions  simultaneously.  The 
.appl  Icatlon  of  the  c-amera  to  the  study  of  flexural  strains  in  bent 
plates  is  shown,  and  the  dctennin.ation  of  the  complete  state  of  two- 
dimensional  strains  is  also  considered.  The  multiple  image-shearing 
camera  uses  an  interference  phenomena,  but  is  less  demanding  chan 
holographic  interferometry  with  respect  to  vibration  isolation  and  the 
coherence  of  the  light  source.  It  is  superior  to  ocher  speckle 
techniques  in  Chat  the  obtained  fringes  are  of  much  better  quality. 
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49. 


A.  J.  Durelll  and  K.  Rajalah,  "Quasi-square  Hole  With  Optimum  Shape 
In  an  Infinite  Plate  Subjected  to  In-plane  Loading"— January  1979. 

This  paper  deals  with  the  optimization  of  the  shape  of  the  comers 
and  sides  of  a  square  hole,  located  In  a  large  plate  and  subjected 
Co  In-plane  loads.  Appreciable  disagreement  has  been  fo\ind  between 
the  results  obtained  previously  by  other  Investigators.  Using  an 
optimization  technique,  the  authors  have  developed  a  quasi-square 
shape  which  Introduces  a  stress  concentration  of  only  2.54  In  a 
uniaxial  field,  the  comparable  value  for  the  circular  hole  being  3. 

The  efficiency  factor  of  the  proposed  optimum  shape  Is  0.90.  whereas 
the  one  of  the  best  shape  developed  previously  was  0.71.  The  shape 
also  Is  developed  that  minimizes  the  stress  concentration  In  the 
case  of  biaxial  loading  when  the  ratio  of  blaxlallty  Is  1:-1. 

50.  A.  J.  Ourelll  and  K.  Rajalah,  "Optimum  Hole  Shapes  In  Finite  Plates 
Under  Uniaxial  Load ,  "—February  1979. 

This  paper  presents  optimized  hole  shapes  in  plates  of  finite  width 
subjected  to  uniaxial  load  for  a  large  range  of  hole  to  plate  widths 
(D/V)  ratios.  The  stress  concentration  factor  for  the  optimized 
holes  decreased  by  as  much  as  44Z  when  compared  to  circular  holes. 
Simultaneously,  the  area  covered  by  the  optimized  hole  Increased 
by  as  much  as  26Z  compared  to  the  circular  hole.  Coefficients  of 
efficiency  between  0.91  and  0.96  are  achieved.  The  geometries  of 
the  optimized  holes  for  the  D/V  ratios  considered  are  presented  in 
a  form  suitable  for  use  by  designers.  It  is  also  suggested  that 
the  developed  geometries  may  be  applicable  to  cases  of  rectangular 
holes  and  to  the  tip  of  a  crack.  This  Information  may  be  of 
Interest  In  fracture  mechanics. 

51.  A.  J.  Durelll  and  K.  Rajalah,  "Determination  of  Strains  in 
Photoelastic  Coatings,”' — Kay  1979 

Photoelastic  coatings  can  be  cemented  directly  to  actual  structural 
components  and  tested  under  field  conditions.  This  Important  advantage 
has  made  them  relatively  popular  In  industry.  The  information  obtained, 
however,  may  be  misinterpreted  and  lead  to  serious  errors.  A  correct 
Interpretation  requires  the  separation  of  the  principal  strains  and  so 
far,  this  operation  has  been  found  very  difficult.  Following  a  previous 
paper  by  one  of  the  authors,  it  is  proposed  to  drill  small  holes  in  the 
coating  and  record  the  birefringence  at  points  removed  from  the  edge  of 
the  holes.  The  theoretical  background  of  the  method  is  reviewed;  the 
technique  necessary  to  use  it  is  explained  and  two  applications  are 
described.  The  precision  of  the  method  is  evaluated  and  found  satisfactory 
in  contradiction  to  Information  previously  published  In  the  literature. 
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52. 


A.  J.  Durelli  and  K.  Rajalah,  "Optimized  Inner  Boundary  Shapes  in 
Circular  Rings  Under  Diametral  Compression, "—June  1979. 

Using  a  method  developed  by  the  authors,  the  configuration  oi  the  inside 
boundary  of  circular  rings,  subjected  to  diametral  compression,  has  been 
optimized,  keeping  cleared  the  space  enclosed  by  the  original  circular 
inside  boundary.  The  range  of  iiameters  studied  was  0.33  <  ID/OD  <0.7. 

In  comparison  with  circular  rings  of  the  same  ID/OD,  the  stress  concentra¬ 
tions  have  been  reduced  by  about  30%,  the  weight  has  been  reduced  by  about 
10%  and  coefficients  of  efficiency  of  about  0.96  have  been  attained.  The 
maximum  values  of  compressive  and  tensile  stresses  on  the  edge  of  the  hole, 
are  approximately  equal,  there  are  practically  no  gradients  of  stress 
along  the  edge  of  the  hole,  and  sharp  comers  exhibit  zero  stress.  The 
geometries  for  each  ID/OD  design  are  given  in  detail. 

53.  A.  J.  Durelli  and  K.  Rajaiah,  "Lighter  and  Stronger," — February  1980. 

A  new  method  has  been  developed  that  permits  the  direction  design  of  shapes 
of  two-dimensional  structures  and  structural  components,  loaded  in  their 
plane,  within  specified  design  constrains  and  exhibiting  optimum  distribution 
of  stresses.  The  method  uses  photoelasticity  and  requires  a  large  field 
diffused  light  polariscope.  Several  problems  of  optimization  related  to 
the  presence  of  holes  in  finite  and  infinite  plates,  subjected  to  uniaxial 
and  biaxial  loadings,  are  solved  parametrically. 

Some  unexpected  results  have  been  found:  1)  the  optimum  shape  of  a  large 
hole  in  a  bar  of  finite  width,  subjected  to  uniaxial  load,  is  "quasi"  square, 
but  the  transverse  boundary  has  the  configuration  of  a  "hat";  2)  for  the 
small  hole  in  the  large  plate,  a  "barrel"  shape  has  a  lower  s.c.f.  than 
the  circular  hole  and  appreciably  higher  coefficient  of  efficiency;  3)  the 
optimum  shape  of  a  tube,  subjected  to  diametral  compression,  has  small 
"hinges"  and  is  much  lighter  and  stronger  than  the  circular  tube.  Applications 
are  also  shown  to  the  design  of  dove-tails  and  slots  in  turbine  blades  and 
rotors,  and  to  the  design  of  star-shaped  solid  propellant  grains  for  rockets. 

54.  C.  Bremond  and  A.  J.  Durelli,  "Experimental  Analysis  of  Displaceir;>.-nts  and 
Shears  at  the  Surface  o-f  Contact  Between  two  Loaded  Bodies," — July  1980. 

The  dlsplacement.s  which  exist  at  the  contact  between  tw’o  loaded  bodies  depend 
on  the  geometry  of  the  surface  of  contact,  the  type  of  the  loading  and  the 
property  of  the  materials.  A  method  has  been  developed  to  deterr'.ine  these 
displacements  experimentally.  A  grid  has  been  photographically  printed  on 
an  interior  plane  of  a  transparent  model  of  low  modulus  of  elasticity.  Th.e 
displacements  were  recorded  photographically  and  the  analysis  was  conducted  on 
the  photographs  of  the  deformed  grids.  Shears  were  determined  from  the  change 
in  angles.  The  precision  of  the  measurements  at  the  Interface  is  estimated 
to  be  plus  or  minus  O.Obnun.  Examples  of  application  are  given  for  the  cases 
of  loads  applied  normally  and  tangentially  to  a  rigid  cylindrical  punch 
resting  on  a  semi-infinite  solft  plate.  Important  observations  can  be  made 
on  the  zones  of  friction  and  of  slip.  The  proposed  method  is  three-dimensional 
and  the  distributions  can  be  obtained  at  several  Interior  planes  by  changing 
the  position  of  the  plane  of  the  grid.  The  limitations  of  the  method  arc 
pointed  out.  The  possibility  of  using  gratings  (12  to  40  Ipmm)  is  considered 
as  well  as  the  advantages  of  using  moire  to  analyze  the  displacements. 
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M.  Erickson  uiui  .  J.  Ihirelli,  "Sf  r^■f  s  i'l  ;  t  r1hui on  Arcu.nd  a  Circular 
Hole  in  Square  I'l  ites,  Ii  nieil  L’niTorir.]  •  in  i  !,<  Cl.uu,  oi,  two  Opi>csiro 
Sides  of  the  .'Square" — Tin  i'i'P'dI  .  t  u  st.ris.-  c  r c  r  i  lu  t  iou  around  ■  (Tr'.ilar 
hole,  located  in  the  cenler  of  square  plaU',  ha.-:  been  deterrri;  c, 
photoela.st  leal  1  V  fpr  llu-  e.ase  of  the  plate  loadid  unlfomly  on  two 
o;-;o.nit<  side.'-.  Tlie  sturly  was  conducted  parar  f- 1  r  i  ca  1 1  y  for  a  larpe  rar.ge 
of  tin  r.ilio  o:  tin'  side  of  ihe  ;;qu.irc  to  the  dian-i  ter  of  thr  Ic  It'.  The 
resuli.s  Lbtaincd  peiirit  the  deteririiual  ion  of  tb.e  stresses  for  a.v.'  biaxial 
condition  and:  verify  a  previ('us  so at  ion  fbtnini'd  for  the  cast  c‘.  the 
pressuriaeJ  hoi  ■.  'I’in  t-;-'.per  i  r'l-nt  al  p  rocedtiT  e  i  .s  br'efly  dc.ser  ibed  . 

56.  A.  J.  Durelll,  M.  Erickson  and  K.  Rajaiah,  "  "This  paper  presents  the 

shapes  that  will  optimize  the  stress  distribution  about  central  holes  in 
square  plates  subjected  to  uniform  load  on  two  opposite  sides  of  the  plate. 

The  study  is  conducted  for  a  large  range  of  hole  to  plate  widths  ratios  (D/W) . 
The  stress  concentration  factor  for  the  optimized  holes  decreased  by  as 
much  as  21%  when  compared  to  the  one  associated  with  a  circular  hole. 
Simultaneously,  the  weight  of  the  plate  with  optimized  hole  is  reduced  by 
as  much  as  36%  as  compared  to  the  circular  hole.  Coefficient  of  efflciecy 
of  around  0.92  is  achieved  for  all  D/W  ratios.  The  geometry  of  the 
optimized  holes  are  presented  in  a  form  suitable  for  use  by  designers. 
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Stress  Distribution  Around  a 
Circular  Hole  in  Square  Plates, 
Loaded  Uniformly  in  the  Plane, 
on  TVo  Opposite  Sides  of  the 
Square 


M.  Erickson^  and  A.  J.  DurallP 


The  complete  stress  distribution  around  a  circular  hole,  located  in 
the  center  of  a  square  plate,  has  been  determined  photoelastically 
for  the  case  of  the  plate  loaded  uniformly  on  two  opposite  sides  The 
study  was  conducted  parametrically  for  a  range  of  the  ratio  D/W  of 
the  diameter  of  the  hole  to  the  side  of  the  square  from  0.20  to  0.S3 
The  results  obtained  permit  the  determination  of  the  stresses  for 
any  biaxial  cortdition  and  verify  a  previous  solution  obtained  for  the 
case  of  the  pressurized  hole.  The  experimental  procedure  is  briefly 
described. 

Introduction 

The  cleuical  problem  of  the  ttrese  distribution  around  a  circular 
hole  in  an  infinite  plate  subjected  to  a  uniaxial  uniform  loading  in  the 
plane  of  the  plate  was  solved  by  Kirsch  [1|  in  closed  form.  The  ap¬ 
preciably  more  complicated  case  of  the  finite  plate  with  the  circular 
hole  was  solved  by  Howland  [2|  using  an  infiniu  series  solution,  but 
results  were  evaluated  only  for  D/W  <  0.5,  D  being  the  diameter  of 
the  hole  and  W  the  width  of  the  plate.  The  distribution  of  stress  fur 


the  cases  when  D/W  >  0.5  was  obtained  experimentally  by  Wahl  and 
Beeuwkes  (S).  The  stress-concentration  factors  referred  to  both  the 
gross  area  and  the  net  area,  for  the  total  range  of  D/W  values  are  given 
in  [4].  The  case  of  the  very  large  hole  in  the  plate,  when  D/H  ap¬ 
proaches  one  presented  some  problems  of  interpretation,  which  have 
been  dealt  with  in  |5]. 

The  stress  distribution  for  the  case  of  a  square  plate  with  a  circular 
hole  was  solved  experimenully  [e]  when  a  uniform  pressure  is  applied 
inside  the  hole,  or  what  is  equivalent  (7|,  when  the  four  sides  of  the 
square  are  subjected  to  uniform  pressure.  The  problem  of  the  square 
plate  with  a  circular  hole,  subjected  to  in-plane  uniform  pressure 
applied  to  two  opposite  sides  of  the  plate,  has  not  been  solved.  The 
problem  is  important  and  if  the  solution  were  available,  the  soluthm 
of  the  previously  mentioned  problem  for  any  ratio  of  biaxiality  could  ^ 
be  obtained  as  a  special  case  by  superposition.  That  is  the  problem 
dealt  with  in  this  Note.  The  solution  is  obtained  photoelastically  for 
a  range  of  D/W  values  from  0.20  to  0.8;1. 

Test  Procedure 

The  analysis  was  conducted  in  a  3-in-sq.  i-in-thick  (Homaliie  100) 
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Fig.  I  Typical  laechremaHc  pallam  around  a  dreular  hola  In  a  aquara  plala 
tubjadad  lo  unNocm  proaaura  on  two  oppoaHa  aMot 
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Fig.  4  nrMS  dMributlom  on  oulor  boundary  o4  a  aquara  plala  wHh  a  round 
hola  aubladad  lo  unMorm  pcooauta  appHad  to  Iwo  oppooHa  oMoa  ol  Iho  plala 
(avaraga  oiraoa  on  Hia  nal  aacllon  uaad  lor  camparim) 


apecimen  (Fig.  1 ).  The  uniform  pressure  is  applied  by  means  of  a 
specially  built  device  as  described  in  |8|.  Two  rubber  hoses,  one  placed 
on  each  of  the  opposite  sides  of  the  plate  are  used.  The  deformation 
of  the  pressurized  hose  is  restrained  by  Plexiglass  sheets.  The  loading 
frame  had  to  be  calibrated  to  determine  the  amount  of  pressure  ac¬ 
tually  applied  to  the  specimen.  For  this  purpose,  a  strain  gaged  load 
cell  was  specially  designed.  The  average  fringe  order  was  computed 
using  the  applied  pressure  and  the  fringe  value  of  the  material  and 
a  check  obtained  by  algebraically  averaging  the  areas  above  and  below 
the  zero  axis  for  those  specimens  with  high  D/M'. 

Seven  specimens  were  used  with  the  inner  hole  diameter  varying 
from  0.6  in.  to  2.5  in.  giving  D/W  values  from  0.2  to  0.83,  where  D  is 
the  hole  diameter  and  tV  the  width  of  the  specimen.  Dark  field  and 
light  field  photographs  were  taken  in  a  diffused  light  polariscope  of 
the  seven  specimens,  subjected  to  pressure  sufficient  to  produce  a 
maximum  of  about  5  fringes  (Fig.  2).  Fractional  fringe  orders  were 
recorded  using  Tardy's  method  of  compensation  at  every  10*  at  the 
edge  of  the  hole  from  0*  (horizontal)  to  90°  (vertical).  Readings  were 
also  taken  on  the  outer  edge  of  the  plate  at  the  0°  and  90*  points.  A 
calibration  test  on  a  2.5-in-dia  round  disk  of  the  material  gave  a  ma¬ 
terial  constant  of  156  Ib/in./fringe. 

The  results  obtained  are  given  as  stress  distributions  along  the 
inside  and  the  outside  boundaries  (Figs.  3  and  4),  and  as  stress  con¬ 
centrations  at  the  intersection  of  the  longitudinal  and  transverse  axes 
with  the  boundary  of  the  hole.  All  values  are  given  parametrically  as 
functions  of  D/W.  These  results  permit,  by  superposition,  the  de¬ 
termination  of  stresses  for  any  ratio  of  biaxial  loading  of  the  plate.  The 
case  of  equal  biaxiality  was  computed  and  is  shown  in  Fig.  5.  The 
values  obtained  verify  those  previously  published  for  the  case  of  the 
hydrostatically  loaded  hole  |6],  using  the  transformation  explained 
in  17]. 

It  may  be  noted  that,  in  the  present  problem,  Ki  increases  as  D/W 
increases  while  for  the  case  of  circular  holes  in  long  rectangular  plates, 
Ki  decreases  as  D/W  increases. 
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Abstnct— This  paper  presents  the  shapes  that  will  optimize  the  stress  distribution  about  central  holes  in 
square  plates  subjected  to  uniform  load  on  two  opposite  sides  of  the  plate.  The  study  is  conducted  for  a 
large  range  of  hole  to  plate  widths  ratios  [DIW).  The  stress  concentration  factor  for  the  optimized  holes 
decreased  by  as  much  as  21%  when  compared  to  the  one  associated  with  a  circular  hole.  Simultaneously, 
the  weight  of  the  plate  with  optimized  hole  is  reduced  by  as  much  as  36%  as  compared  to  the  circular  hole. 
Coefficient  of  efficiency  of  around  0.92  is  achieved  for  ii\  DIW  ratios.  The  geometry  of  the  optimized  holes 
are  presented  in  a  form  suitable  for  use  by  designers. 

INTRODUCTION 

This  paper  is  one  in  a  series  of  papers  dealing  with  the  optimization  of  discontinuities  in  two 
dimensional  stress  fields.  The  optimum  shape  of  a  hole  in  an  infinite  plate  subjected  to  uniaxial 
uniformly  distributed  load  was  presented  in[l].  The  shapes  to  be  given  to  an  optimized  central 
hole  in  finite  plates  subjected  to  uniaxial  uniform  loading  have  been  given  in  [2]  for  different 
ratios  of  the  diameter  of  the  hole  to  the  width  of  the  plate.  The  optimized  inner  boundary 
shapes  of  rings  with  circular  outer  boundaries  subjected  to  diametral  compression  have  been 
given  in  [3]  for  different  ratios  of  outer  to  inner  diameters.  The  description  of  the  basic  features 
of  the  method  have  been  presented  inf4, 5],  References  to  other  contributions  in  the  literature 
can  be  found  in  [5]  among  which  the  most  important  one  is  due  to  Hey  wood  [6]. 

OPTIMIZATION  PROCEDURE 

The  method  consists  in  using  photoelasticity,  in  a  systematic  way.  to  idealize  a  configuration 
so  that  its  boundaries  do  not  have  gradients  of  stress  along  the  length  of  the  boundaries.  In 
other  words,  the  structure  will  have  stresses  uniformly,  or  almost  uniformly  distributed  along 
the  boundary.  The  procedure  permits  the  direct  design  of  the  geometry  of  the  structure  rather 
than  the  conventional  step  by  step  design  and  analysis,  satisfying  the  requirement  that  the 
maximum  stress  should  be  lower  than  an  allowable  stress  and  at  the  same  time,  the  distribution 
should  be  as  efficient  as  possible.  The  geometric  constrains  for  the  problem  are  stripulated 
initially.  A  transparent  model  of  the  structure  is  placed  in  a  diffused  light  circular  polariscope. 
(The  material  of  the  model  should  exhibit  birefringence  when  under  load  and  should  be 
sufficiently  sensitive  to  produce  several  fringes  of  interference.)  The  operator  should  be  able  to 
work  on  the  model  with  a  hand  file  or  portable  router  at  the  same  time  that  he  looks  at  the 
model  through  the  analyzer.  The  method  requires  that  the  operator  file  away  material  from  the 
boundaries  starting  at  the  points  where  the  stress  (and  therefore  the  fringe  order)  is  at  a 
minimum.  The  filing  operation  redistributes  the  fringes.  The  operator  continues  filing  away 
material  from  the  low  stressed  zones  of  the  boundary  until  as  much  as  possible  of  the  length  of 
the  boundary  shows  the  same  order  of  birefringence.  This  is  easy  to  detect  because  at  that 
moment,  the  same  fringe  falls  all  along  the  length  of  the  boundary.  If  the  body  has  more  than 
one  boundary,  it  may  be  necessary  to  operate  by  steps  back  and  forth  from  one  boundary  to  the 
other.  In  some  cases,  white  light  may  be  more  practical  than  monochromatic  light,  the  objective 
being  then  to  have  the  same  color  along  the  boundary. 

The  degree  of  optimization  can  be  evaluated  quantitatively  by  a  coefficient  of  efficiency; 

fSi  ^  r*! 

,  \  <r,  ds  \  a,  ds 
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where  a,  is  the  tangential  stress,  represents  the  maximum  allowable  stress  (the  positive  and 
negative  superscripts  referring  to  tensile  and  compressive  stresses,  respectively).  So  and  5|  are 
the  limiting  points  of  the  segment  of  boundary  subjected  to  tensile  stresses  and  Si  and  S2  are 
the  limiting  points  of  the  segment  of  boundary  with  compressive  stresses.  The  significance  of 
the  coefficient  of  efficiency  was  discussed  in  Refs.  [2,  S].  The  above  criterion  will  be  used  in  the 
present  work  to  evaluate  the  optimized  hole  shapes.  The  design  procedure  will  be  particularly 
useful  for  components  made  with  brittle  materials,  or  components  made  with  ductile  materials 
subjected  to  fatigue. 
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Pig.  2.  Typical  isochromatic  pattern  around  an  optimized  hole  in  a  square  plate  subjected  to  a  uniaxial 

uniform  pressure. 
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Optimum  shapes  of  central  holes  in  square  plates 

In  this  paper,  using  the  technique  described  above,  optimized  hole  shapes  are  presented  in 
square  plates  subjected  to  uniaxial  uniform  compression. 

EXPERIMENTAL  DETAILS 

Experiments  were  conducted  with  0.23  in.  (5.8  mm)  thick  Homalite-100  plates  (fringe  con¬ 
stant  of  1561b/in-fr  (27.0kN/m-fr)).  The  plate  size  was  chosen  as  3  x3in.  (76.2  x  76.2  mm)  for 
all  DIW  ratios.  Optimization  was  carried  out  for  OW  =  0.20;  0.41;  0.56;  0.65;  0.71;  0.79  and 
0.89,  with  the  models  subjected  to  uniaxial  uniformly  distributed  compression.  Material  was 
removed  from  the  low  stress  regions  by  careful  hand  filing.  To  improve  the  precision,  in 
particular  at  the  corner  zones,  the  operator  used  a  binocular  magnifier  with  a  set  of  polarizer 
and  quarter  wave  plates  attached  to  each  of  its  lenses,  during  the  filing  process.  The  uniform 
compression  on  the  two  parallel  boundaries  was  applied  following  the  methods  developed 
previously  [7, 8]  and  used  recently  in  (9).  The  position  of  a  pressurized  rubber  tube  located  in 
the  special  device  is  shown  in  Fig.  I.  This  loading  frame  had  to  be  calibrated  to  determine  the 


Fig.  3.  Sirens  ratios  for  points  on  the  boundary  of  optimized  circular  holes  in  a  square  plate  subjected  to  uniaxial 

uniform  pressure. 


Fig.  4.  Comparison  of  stress  concentration  factors  for  an  optimized  hole  vs  a  circular  hole  in  a  square  plate 
subjected  to  uniaxial  uniform  pressure. 
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amount  of  pressure  actually  applied  to  the  specimen.  For  this  purpose,  a  strain  gaged  load  cell 
was  specially  designed. 


RESULTS 

The  isochromatic  pattern  for  a  typical  hole  shape  is  shown  in  Fig.  2.  The  stress  distributions 
around  circular  and  optimized  holes  for  the  £VW  ratios  considered  are  presented  in  Fig.  3.  The 
s.c.f.  for  the  tensile  and  compressive  regions  of  the  circular  and  optimized  holes  for  different 
D/W  ratios  are  plotted  in  Fig.  4.  The  information  on  coefficient  of  efficiency  and  percentage 
weight  reduction  achieved  are  given  in  Fig.  5.  The  stress  distributions  on  the  outer  edges  are 
shown  in  Fig.  6. 


Fig.  S.  Percentage  weight  reduction  and  coefficient  of  efficiency  for  an  optimized  hole  in  a  square  plate 
subjected  to  a  uniaxial  uniform  pressure. 


Fig.  6.  Stress  distribution  along  the  outer  boundaries  of  a  square  plate,  with  an  optimized  hole,  subjected  to 
uniform  pressure  on  two  opposite  sides. 
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Fi(t.  Radii  of  the  elements  of  holes  producing  optimum  distriholion  of  stress  in  square  plates  subjected  to 

uniaxial  uniform  pressure. 


The  experimentally  developed  optimum  hole  geometries  have  been  fitted  w  ith  a  combination 
of  circles  of  different  diameters  and  common  tangents  at  the  points  of  intersection.  The  hole 
geometries  for  the  different  DIW  ratios  are  given  in  Fig.  7.  The  information  given  in  Fig.  7  is 
consolidated  in  Fig.  8  in  a  graphical  form  to  make  easy  its  application. 

Dl.SCtlSSION 

The  isochromatic  pattern  in  Fig.  2  shows  that  the  newly  developed  hole  shape  exhibits  a 
high  degree  of  optimization  with  the  stresses  remaining  uniform  along  major  portions  of  the 
tensile  as  well  as  compressive  segments  of  the  boundary.  The  information  in  Figs.  }  and  4 
shows  that,  as  compared  to  the  circular  holes,  the  optimum  shapes  have  lead  to  significant 
reduction  in  s.c.f.,  the  reduction  ranging  from  16%  for  DIW  =  0.14  to  about  21%  for  DIW  - 
0.84. 

A  coefficient  of  efficiency  of  about  0.92  has  been  achieved  for  all  D/W  ratios,  as  seen  from 
Fig.  5.  The  optimum  shapes  have  also  lead  to  significant  reduction  in  weight  as  compared  to  the 
circular  holes,  in  the  range  of  1%  for  DIW  =  0.14  to  36%  for  DIW  =  0.84. 

Comparison  of  the  optimized  shapes  and  the  s.c.f.  presented  here  for  square  plates  with 
tho.se  given  in  Ref.  (2)  for  long  rectangular  plates  under  identical  loading  conditions  (except  for 
the  reversal  of  sign)  brings  out  certain  interesting  features. 

For  DIW  <0.3,  the  optimized  shapes  are  found  to  be  identical,  whereas  for  DIW  ''  0.3.  the 
two  shapes  tend  to  become  quite  different  especially  on  the  hole  edges  perpendicular  to  the  load 
axis.  As  DIW  becomes  larger,  the  hole  shapes  in  square  plates  are  predominantly  influenced  b\ 
the  bending  effect  of  the  horizontal  segments.  No  such  bending  effect  can  be  observed  in  long 
rectangular  plates. 

For  a  very  large  value  of  DIW,  the  problem  can  be  considered  as  that  of  a  portal  frame 
subjected  to  uniform  load  on  top.  The  optimum  shape  presented  can  be  taken  as  the  shape  to  be 
given  to  a  stress  optimized  portal  frame. 

While  the  s.c.f.  for  square  plates  decreases  initially  slightly  and  then  increases  monotonic- 
ally  with  increase  in  DIW  (Fig.  4).  it  was  found  to  decrease  monotically  for  long  rectangular 
plates  as  DIW  increases.  For  a  given  DIW  ratio  beyond  DIW  =  0.5.  the  s.c.f.  for  square  plates 
is  found  to  be  very  much  higher  than  that  for  long  rectangular  plates.  The  reason  for  this 
increase  can  be  traced  to  the  behavior  of  the  square  plate  as  a  portal  frame,  which  increases 
appreciably  the  bending  stresses  in  both  the  vertical  and  horizontal  members. 

For  large  DIW.  (he  vertical  elements  of  the  plate  become  thinner  and  the  fringes  are  not 
only  parallel  to  the  inner,  but  also  to  a  large  extent,  to  the  outer  boundaries  indicating  a  linear 
variation  of  stress  across  the  section,  li  also  shows  (hat.  for  large  DIW.  as  the  inner  edge  gets 
optimized,  the  outer  edge  also  tends  to  become  optimum.  The  distribution  of  .tresses  on  the 
outer  edges  shown  in  Fig.  6  confirms  this.  It  is  also  seen  from  Fig.  6  that  the  stress  distribution 
becomes  favorable  also  on  the  loaded  edges  with  the  tangential  stresses  remaining  constant 
over  a  considerable  length. 

It  is  also  seen  from  Fig.  4  that  the  location  of  maximum  stress  shifts  from  the  vertical  edge 
to  the  horizontal  edge  for  CVW>0.81. 
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Optimization  of  Hole  Shapes  in  Circular  Cylindrical 
Shells  under  Axial  Tension 


Optimization  of  hole  shapes  in  cylindrical  shells  under 

axial  tension  is  obtained  using  two-dimensional  photoelastic 

and  birefringent-coating  techniques,  resulting  in  uniform  membrane  stresses 

along  the  tensile  and  compressive  segments  of  the  hole  boundary  ' 

by  K.  Rajaiah  and  A.J.  Durelli 


ABSTRACT — Hole  snapes  are  opumized  m  circulai  cylindrical 
sneiis  suPjecied  to  a«ia'  load  considering  oniy  Ihe  pieoominanily 
large  membrane  stresses  pieseni  around  me  holes  Two- 
dimensional  pholoeiastic  isochiomatics  obtained  with  a 
special-purpose  polanscope  are  utilized  lor  the  optimization 
process  The  process  leads  Ip  a  sigmlicanl  decrease  in  me 
membrane  stress-concentration  laclor  and  a  modest  decrease 
in  weight,  thus  yielding  a  considerable  increase  in  sirength-io- 
weight  ratio  This  paper  presents  results  lor  certain  typical 
ratios  ol  hole  diameter  to  shell  diameter  Previous  theoretical 
and  experimental  studies  lor  the  circular  hole  have  also  been 
verilied 

Introduction 

Optimization  of  hole  shapes  in  engineering  structures 
like  plates  and  shells  is  important  from  fatigue  and 
minimum-weight  considerations.  In  a  series  of  papers, 
Durelli  and  his  associates'"*  have  recently  presented  a 
practical  way  of  arriving  at  optimum  hole  shapes  in  uni- 
axially  loaded  plates  and  rings  from  simple  two-dimen¬ 
sional  photoelastic  experiments  by  removal  of  material 
from  low-stress  regions  around  the  hole  and  making  an 
isochromatic  fringe  coincide  separately  with  the  tensile 
and  compressive  segments  of  the  boundary.  This  process 
leads  to  a  significant  reduction  in  the  stress-concentration 
factor  as  well.  The  present  paper  is  an  extension  of  the 
same  concepts  to  the  optimization  of  hole  shapes  in 
circular  cylindrical  shells  under  uniaxial  tension. 

The  stress  pattern  around  holes  in  cylindrical  shells  is 
complicated  by  the  fact  that  both  membrane  and  bending 
stresses  are  present.  Fortunately,  for  circular  holes  in 
circular  cylindrical  shells  under  axial  tension,  Peterson's’ 
handbook  shows  that  membrane  stresses  (computed  by 
Van  Dyke)  are  predominant  on  the  hole  edge  while 
bending  stresses  are  of  a  much  smaller  magnitude  (Fig.  I). 
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As  a  percentage  of  the  total  stress,  the  bending  stresses 
constitute  less  than  16  percent  for  t/R  >  0.004  (Fig.  2). 


Fig.  1  — Membrane  and  bending  stress-concentration 
factors  (or  circular  holes  in  circular  cylindrical  shell 
under  axial  tension 
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Fig  2— Percentage  ol  maximum  bending  stress  present 
around  the  hole  lor  various  curvatures  d  ol  the  hole 
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In  an  ideal  siiuaiion,  one  would  aiiempt  to  opiimi;e  ihc 
hole  shape  considering  both  the  membrane  and  ihc 
bending  stresses.  Obviously,  such  an  exercise  would  be 
quite  involved.  From  a  practical-design  point  of  view,  it 
should  be  enough  if  the  holt  shapes  are  optimieed  by 
taking  into  account  only  the  predominantly  large  mem¬ 
brane  stresses.  The  method  proposed  in  this  paper  is 
based  on  this  idea. 

Defining  an  optimised  hole  as  one  in  which  uniiorm 
membrane  stresses  are  present  in  both  the  tensile  and 
compressive  segments  of  the  hole  boundary,  it  is  found 
that  two-dimensional  photoelast icily  (either  by  trans¬ 
mission  or  by  reflection)  can  be  used  to  obtain  optimized 
hole  shapes  using  shell  models  made  of  birefringent 
materials.  In  the  present  investigation,  this  technique  is 
utilized  for  the  optimization  of  hole  shapes  in  circular 
cylindrical  shells  under  uniaxial  loading.  Optimized  hole 
shapes,  the  corresponding  stress-concentration  factor  and 
stress  distributions  around  the  holes  are  presented  for 
some  typical  ratios  of  hole  diameter/shell  diameter. 
The  geometric  constrain  specified  for  the  optimized  hole 
is  that  the  optimized  boundary  should  be  in-between  the 


Fig  3— Loading  fixture 
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Fig.  4— Transmission-polariscope  arrangement 


circle  of  diameter  2o,  and  the  square  of  side  2a. 

Model  Making 

Three  different  types  of  shells  were  made  Preliminarv 
experiments  used  commercially  available  Perspex  tubes  ol 
3-in.  (76-mm)  nominal  diameter  with  0. 125-in.  (jl-mmi 
nominal  wall  thickness  (identified  as  Type  I  shell).  Vkhilc 
the  thickness  distribution  was  within  S  percent  of  the 
nominal  value,  the  cylinder  was  not  sufficiently  straight. 
Due  to  this,  certain  bending  effects  were  observed  The 
advantage  of  this  material  was  the  negligible  time-edge 
stresses  around  the  hole,  but  the  major  problem  was  the 
large  load  required  to  produce  a  sufficient  number  ol 
fringes. 

The  second  set  of  expieriments  used  shells  made  of 
epoxy  sheets,  1.3  mm  ±0.0S  mm  thick  cast  between  two 
substantially  thick  flat  Perspex  sheets  (identtfied  as  Type 
II  shell).  Each  sheet  in  a  semipolymerized  state  was  cut  to 
size  and  wrapped  around  a  SCi-mm-diani  ground-metal 
cylinder  coaled  with  a  releasing  agent.  The  edges  of  the 
epoxy  sheet  touching  each  other  were  joined  together  with 
a  5-mm-wide,  very-thin  strip  providing  a  bonded  butt  joint. 
To  provide  symmetry  to  the  shell,  an  identical  strip  was 
bonded  on  the  diagonally  opposite  side.  The  open  ends 
of  the  cylindrical  shell  were  reinforced  over  a  distance  of 
1.0  in.  (25  mm)  with  bonded  epoxy  strips  to  provide 
stiffness  near  the  end  fixture.  After  complete  curing,  the 
cylinder  was  removed  and  tested.  For  a  no-hole  condition, 
it  was  verified  that  there  were  no  residual  stresses  except 
very  close  to  the  joints,  and  that  the  reinforcement  affects 
the  stress  distribution  over  the  uniform  section  only  to  a 
distance  of  about  twice  the  shell  thickness.  Even  with  the 
introduction  of  holes,  there  were  no  measurable  effects 
of  the  reinforcing  strips  over  the  hole  stress  distribution. 
The  advantage  of  this  technique  was  the  fine  transmission 
properties  of  the  shell  and  absence  of  machining  process. 

The  last  set  of  experiments  were  carried  out  with 
integrally  cast  circular  cylindrical  shells  (identified  as 
Type  111  shell)  using  a  plaster -of-paris  mold  with  a  centrally 
located,  slightly  tapered  metal  cylinder  coated  with  a 
silicone-grease-releasing  agent.  After  complete  curing,  the 
inner  metal  cylinder  was  first  pushed  out  using  a  hydraulic 
press  and  then  the  outer  plaster-of-paris  shell  was  removed. 
Both  the  inner  and  outer  surfaces  of  the  epoxy  shell* 
were  turned  on  a  lathe  to  the  required  dimensions  (50-mm 
diameter  and  0.4  mm  thick).  The  surfaces  needed  a  thin 
coating  of  oil  to  make  them  transparent.  Dimensionally, 
these  cylinders  were  the  best.  In  the  last  two  experimenis, 
in  order  to  overcome  the  time-edge  effect  around  the 
holes,  the  experiments  were  carried  out  immediately  after 
the  introduction  of  the  holes.  Optimization  was  also  com¬ 
pleted  on  the  same  day. 

The  introduction  of  the  circular  hole  was  carried  out  in 
different  stages  starting  from  a  small-diameter  hole  and 
progressively  increasing  it  to  the  required  diameter  and 
then  finishing  with  a  reamer.  The  holes  were  stress  free 
for  all  practical  purposes.  The  length  of  all  the  shells  used 
was  sufficient  for  the  ends  to  have  negligible  effect  on  the 
hole.  No  creep  effects  were  noticed  in  any  of  the  shells. 

Loading  System 

Uniform  axial  loading  of  both  tension  and  compression 
types  was  applied  on  the  models.  The  compression  load 
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was  given  through  rubber  sheet  and  wooden  blocks  kept 
on  the  model  with  a  12-mm  ball  located  at  the  center. 
The  shell  could  be  axially  loaded  and  the  uniformity  of 
the  stress  at  some  distance  away  from  the  load  checked. 

The  tensile  load  was  applied  through  shear  v^ith  tightly 
gripped  end  rings  and  a  disk  with  a  spherical-ended  rod  at 
the  center  ensuring  the  alignment  of  the  shell  (Fig  3).  In 
these  tests,  compression  loading  was  more  difficult  to 
control  than  tensile  loading  and  all  the  reported  results 
correspond  to  tensile  loading. 

Type  of  Polariscope  Used 

Conventional  transmission  polariscopes  are  clearly 
unsuitable  for  the  study  of  the  present  problem.  A 
polariscope,  based  on  Houghton's  experiments,*  with  the 
polarizer  and  first  quarter-wave  plate  having  their  planes 
parallel  to  the  plane  of  the  hole  and  located  at  the  center 
of  the  shell,  was  used.  A  diffused  130-W  sodium-vapor 
lamp  was  used  as  light  source  on  one  side  of  the  shell 
while  the  second  quarter-wave  plate  and  the  analyzer  were 
located  on  the  other  side.  The  fringes  were  recorded  using 
a  bellows-type  camera  with  a  25-in.  (62.5-cm)  focal-length 
lens  (Fig.  4). 

For  the  reflection  studies,  a  conventional  polariscope  of 
the  type  shown  in  Fig.  5  was  tried.  With  this  polariscope, 
the  shadow  around  the  hole  edge  is  very  pronounced.  The 
results  obtained  are  not  very  reliable.  Due  to  this  difnculty, 
the  reflection  method  was  only  used  in  the  present  study 
as  an  occasional  check.  A  special  reflection  polariscope  of 
the  type  used  by  Slot'  and,  more  recently,  by  Durelli  e/ 
al.‘  would  appear  to  be  ideal  for  this  problem.  With  one 
polariscope  as  with  the  other,  only  membrane  stresses 
are  recorded. 

Optimization  Process 

For  optimization,  the  constraint  was  stipulated  that  the 
boundary  of  the  hole  lie  between  the  circle  of  diameter 
2o  and  the  square  side  2o.  The  process  followed  is  very 
similar  to  the  one  mentioned  in  Refs.  1-4.  It  starts  with 
the  viewing  of  the  isochromaiics  around  the  circular  hole 
in  the  cylinder  loaded  axially  ,  using  the  polariscope  des¬ 
cribed  above.  From  a  study  of  the  isochromaiics,  the 
(ensile  and  compressive  segments  and  the  low  and  high 
stress  regions  around  the  hole  boundary  are  established. 
The  operator  starts  Filing  away  material  from  the  boundary 
at  the  points  where  the  stress  is  minimum,  say.  on  the 
tensile  segment  of  the  boundary  ,  and  with  the  hole  region 
in  view  through  the  analyzer.  This  operation  redi.stributes 
(he  stresses  bringing  down  the  maximum  stress  along  that 
segment.  The  operation  is  continued  until  an  isochromaiic 
fringe  coincides  with  the  tensile  boundary.  The  same 
operation  is  conducted  on  (he  other  segments. 

Results  and  Discussion 

The  isochromaiic  patterns  for  a  typical  circular  and 
corresponding  optimized  holes  are  shown  in  Fig.  6.  The 
stress  distributions  at  the  boundaries  of  the  holes  for  three 
different  ratios  of  hole  diameter  to  shell  diameter  are 
presented  in  Fig.  7.  For  the  sake  of  comparison.  Fig,  7 
also  includes  the  distributions  corresponding  to  (he 
circular  holes  as  determined  in  the  present  experiments. 
The  empirically  developed  geometry  has  been  fitted  with 
a  combination  of  circles  of  different  diameters  and 
common  tangents  at  the  points  of  intersections.  The 
geometry  of  a  typical  optimized  shape  is  shown  in  Fig.  8. 


Jessop  el  a!.'  have  determined  the  stress-concentration 
factors  in  circular  tubes  with  transverse  circular  holes 
under  tension  for  some  typical  values  of  ^  and  I/R ,  using 
(he  stress-freezing  technique.  Comparison  of  their  results 
with  the  present  values,  as  well  as  with  the  analytical 
results  quoted  in  Ref.  5,  shows  close  agreement.  Fig.  9, 
but  are  about  S  percent  lower. 

The  condition  under  which  the  tests  were  conducted 
required  that  the  diameter  of  the  holes  used  be  small.  This 
increased  the  difficulties  of  the  analysis.  It  is  estimated 
that  the  level  of  precision  obtained  is  5  percent.  The 
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Fig.  6— Isochromatic  fringe  patterns 
around  circular  and  optimized  boles  in  a 
circular  cylindrical  shell  under  axial 
load  (a/ft  =  0.24,  f/R  =  0.056) 


Fig  7— Stress  distribution  at  the  boundary  o(  optimized 
and  circular  holes  in  circular  cylindrical  shells  under  axial 
load((/R  «  0.056) 
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Fig-  8— Optimi2ed  geometries  of  holes  for 
minimum  stress-concentration  factor  in  a  circular 
cylindrical  shell  under  axial  load 


Fig.  9— Stress-concentration  factors  for  a  sirgle 
circular  and  optimued  hole  in  circular  cylindrical 
shells  Subjected  to  axial  loading 


optimization  process  is  also  difnculi  on  small  holes.’ 
However,  the  optimized  hole  shapes  and  the  percentage 
reduction  in  stress-concentration  factor  can  be  considered 
to  be  reliable.  For  better  accuracy,  larger  hole  sizes  in 
large-diameter  shells  should  be  used. 

In  Fig.  9,  comparison  of  the  siress-concentration- 
facior  values  for  circular  and  optimized  holes  shows  that 
the  optimization  process  has  brought  down  the  stress- 
concentration  factor  by  as  much  as  25  percent.  The 
reduction  in  stress-concentration  factor  is  larger  for  the 
large  hole  sizes.  Larger  gains  can  be  anticipated  for  j3  >  I. 


The  stress-concentration  factor  increases  with  in.-reasing 
hole  diameter  for  optimized  holes  also.  The  optimum  hole 
shapes  obtained  for  the  cases  investigated  arc  close  to  a 
double  barrel. 

Comparison  of  the  results  obtained  using  shells  with 
(Type  II)  and  without  (Type  111)  reinforcements  indicate 
that,  for  the  range  of  hole  diameters  considered,  the 
stress-concentration-factor  values  in  both  cases  are  in 
close  agreement.  The  optimization  process,  apart  from 
Ihe  reduction  in  stress-concentration  factor  has  also 
resulted  in  a  decrease  in  weight  as  evidenced  by  the 
increase  in  the  hole  area  in  the  optimized  case.  The 
increase  in  hole  area  is  found  to  be  of  the  order  of  S 
percent  to  II  percent.  One  can  conclude  that  the  optimiza¬ 
tion  process  has  thus  led  to  a  significant  increase  in  the 
strength/weight  ratio  of  the  shell. 

As  mentioned  earlier,  the  hole  optimization  has  been 
carried  out  considering  only  the  effect  of  membrane 
stresses.  Whether  this  process  leads  to  an  increase  or 
decrease  of  the  bending  stresses  around  the  hole  is  not 
known  yet.  For  a  complete  understanding  of  the  problem, 
estimation  of  the  bending  stresses  for  both  (he  circular 
and  optimized  shapes  would  be  required,  particularly  for 
very  thin  shells. 

As  a  matter  of  interest,  the  shape  shown  in  Fig  9  for 
the  small  hole  could  be  compared  to  the  one  developed  in 
Ref.  2,  for  (he  case  of  the  large  plate  under  uniaxial 
loading,  and  the  one  developed  in  Ref.  4  for  the  case  of 
the  large  plate  subjected  to  equal  biaxial  loading  of 
opposite  signs.  It  can  be  noted  (hat  the  transverse  boun¬ 
dary  of  (he  optimized  hole  in  (he  case  of  the  plate  under 
uniaxial  load  is  flat,  rather  than  curved  as  in  the  other 
two  cases.  It  is  possible  that  further  work  on  the  shell 
case  will  show  that  an  improvement  in  optimization  can 
be  obtained  by  using  a  flatter  transverse  boundary  for 
the  hole. 
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Hole  shapes  are  optimized  in  circular  cylindrical  shells  subjected  to  axial 
load  considering  only  the  predominantly  large  membrane  stresses  present  arounc 
the  holes.  Two  dimensional  photoelastic  isochromatlcs  obtained  with  a  special 
purpose  polariscope  are  utilized  for  the  optimization  process.  The  process 
leads  to  a  significant  decrease  in  the  membrane  s.c.f.  and  a  modest  decrease 
in  weight,  thus  yielding  a  considerable  increase  in  strength  to  weight  ratio. 
This  paper  presents  results  for  certain  typical  ratios  of  hole  diameter  to 
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circular  hole.  Coefficient  of  efficiency  of  around  0.92  is  achieved 
for  all  D/W  ratios.  The  geometry  of  the  optimized  holes  are  presented 
in  a  form  suitable  for  use  by  designers. 
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